NBI-HE-92-28 
May 1992 



Non-perturbative 2d quantum gravity and hamiltonians 

unbounded from below 



J. Ambj0rn and C.F. Kristjansen 

The Niels Bohr Institute 
Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 



Abstract 

We show how the stochastic stabilization provides both the weak coupling genus 
expansion and a strong coupling expansion of 2d quantum gravity. The double 
scaling limit is described by a hamiltonian which is unbounded from below, but 
which has a discrete spectrum. 
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1 Introduction 



One of the problems of euclidean quantum gravity is that the action is not bounded 
from below. One might think that this is not a serious problem in two dimensions 
where the Einstein-Hilbert action is topological: 

J d 2 ^R = 27TX (1) 

However, in case one attempts a summation over all genera, as is needed in string 
theory and as might be needed in quantum gravity, it creates ambiguity. This 
ambiguity is highlighted by analyzing 2d gravity in the matrix model approach. 
It is conveniently formulated as follows: In the regularized path integral the class 
of manifolds to be included can be chosen to be the piecewise linear orientable 
manifolds, which one can get by gluing together equilateral triangles. The curvature 
assignment can then be done according to Regge calculus which means that the 
curvature is assigned to the vertices and is proportional to the number of triangles 
sharing a given vertex. The gluing of K labeled oriented triangles can be done in 
a systematic way by associating with each oriented link a hermitian N x N matrix 
(f) al3 (see fig.l) 




(f)-y a 

Figure 1: Labels of a triangle, 
choosing an appropriate gaussian measure 

dfj, N ((j)) = c N exp(-]-J2 |0a/3| 2 ) n rfRe ^II dIni ^' ( 2 ) 

1 a,0 a<(3 a<(3 

and taking the action to be proportional to the product of matrices around the 
oriented triangle: 

-4(0) = 0a/30/3 7 7 a = Tr 3 . (3) 

In this way all closed triangulated surfaces of the kind mentioned will be generated 
by the following field theoretical model 

Z(h,N) = J e hA W'^dn N {<t>). (4) 

In fact, by expanding e M ^, the Wick contractions of the term (h/\^N) K A(<fi) K / K\ 
will generate all oriented 2D simplicial manifolds consisting of K triangles. Further 
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the 1/iV-expansion organizes the surfaces according to powers of N. For a surface 
with Euler number x the total power of N is \. 
If we write the continuum euclidean action as 

S(X, G) = \J d 2 ^ - I d 2 ^R (5) 
a comparison with ([]]) leads to the identification of e X//G and N x : 

£ = log N (6) 

Of course the field theory (Q) is not well defined since the action is not bounded 
from below. It should be emphasized that this unboundedness has nothing to do with 
the cubic action. Had we decided to glue together regular squares the action Tr (p 3 
in (|3|) would have been changed to a quartic action Tr0 4 , but it would appear with 
the wrong sign, since each surface should be assigned a positive weight. However, 
each term in the 1/iV-expansion is perfectly well defined and has a finite radius of 
convergence. The difference between the restricted summation provided by each 
term in the 1/iV-expansion and the unrestricted summation over all surfaces is the 
following: For a fixed topology, i.e. fixed Xi t ne number of surfaces which can be 
constructed by gluing together K equilateral triangles grows exponentially, while 
the total number of such surfaces grows factorially, i.e. much faster: 

Mtot(K) > const. K\ 

The exponential bound in (|7|) leads to a critical point which is the same for all \ 
and which allows us to take a scaling limit for /i — >• /i c for each individual x- O ne 

gets 

Z xip ,G) = ^^M x (K)e->"<, (/t _^ 
In (H) \i is related to h in (|) by 

h = (9) 

The so-called double-scaling limit is an attempt to go beyond this expansion, 
and can be viewed as a formal renormalization of the gravitational constant G in 
front of the Einstein-Hilbert action. The key observation is that 

7x = -^X + 2 (10) 

which allows us to write 

G) ~ £c x e*/ G * (11) 
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where the renormalized gravitational constant is defined as 

^ = ^ + ^log(/^-^) (12) 

The remarkable fact that the partition function in the scaling limit is only a function 
of the renormalized gravitational coupling constant was first observed in [III in 
the context of Liouville theory, and in the context of the matrix model in |3], [5], |J. 

The factorial growth ([^D of the total number of piecewise linear surfaces is re- 
flected in a factorial growth of the coefficients c x in (|11]). The series is only an 
asymptotic series and not even Borel summable since all the coefficients are posi- 
tive. This is consistent with the unboundedness of the action (|3|). The origin of the 
unbounded action (|3p is therefore different from the usual problem caused by the 
conformal mode in euclidean gravity, but it is a problem which will be present even 
more severely in higher dimensional gravity. 

Two possible cures of the unboundedness of the action in euclidian quantum 
gravity were suggested long time ago. One method is based on stochastic quantiza- 
tion and uses the bounded Fokker-Planck potential || . The other method involves a 
rotation of the integration contour in the functional integral into the complex plane 
such that the functional integral converges They were both intended to be used 
to cure the problem of the conformal mode, but have also been applied to regularize 
the summation over topology in 2d gravity [§, [|, [L4|, [10], [L2|, [TJj . The outcome is that 



the two methods do not agree JT5 1 . In fact the contour rotation leads to a complex 



partition function and complex correlators, a signal that either the method is no 
good or the theory itself is incomplete. The stochastic method has no obvious flaws 
and in the following will will analyze this method in detail^]. 

2 Stochastic regularization 

The stochastic quantization scheme for a <i-dimensional euclidean field theory with 
an action S[<p] shows that the vacuum expectation value of any operator Q can be 
interpreted as the ground state expectation value in a (d + l)-dimensional quantum 
theory: 

<Q>=^Jd<J) e- s ^ h Q[<J)} =Jd4> tfgfe] Q[<j>]. (13) 

Here 

e -S[<f,]/m 

tt o [0] = 7=- (14) 



^^Here it should also be mentioned that yet another suggestion of a non-perturbative regular- 
ization can be found in the papers by Morris et al. [jT6f . This regularization has no obvious flaws 
either, but does unfortunately not agree with the stochastic method [0. At the moment we do 
not have any convincing arguments in favor of one of the the two methods. 



4 



is the ground state of the theory, determined by the Fokker-Planck hamiltonian 
Hpp = J d d x 

This hamiltonian is a positive semidefinite operator and can be written as 

H FP = J d d x R\x)R{x), R = ^ + ^ (16) 

and it is readily verified that (|T^) is an eigenvector of Hpp corresponding to energy 
E = 0. 

For this to make sense it has been assumed that the action S[<f>] produces a 
normalizable wavefunctional by fll4[). This is clearly not the case if S[<j)} is unbounded 
from below. Formally the wavefunctional (|14]) will still satisfy Hpp^ Q = 0, even for a 
bottomless action, but it does not qualify as a groundstate as it is not normalizable. 
The true ground state will correspond to the lowest eigenvalue of Hpp where the 
eigenstate is normalizable. If we denote these by Eq and (again) by \l/o, we can 
write in the case of a bottomless action: 

-s eff m/2h 

%[4>] = , S eff ^S, ^ T) >0. (17) 

This equation can serve as a definition of a new stabilized action. The expectation 
values of observables are defined as in (Jl3|): 

(Q) = -J— [ d<j> Q[<j>] e- s ^' h = <vfo|Q|vI>o> (18) 

and it can be shown that they have the same classical limit, the same perturbative 
expansion in coupling constants and the same 1/N expansion as in the ill-defined 
bottomless theory. Especially the last point will be important for us. 

Let us now apply the above formalism to the matrix models describing 2d-gravity. 
The partition function of 2d-gravity is in the matrix model approach given by 

Z = Jd^e- NTlV ^^ (19) 

where (f> denotes a. N x N hermitean matrix and the potential V(x) is characterized 
by starting with a quadratic term and then containing a number of wrong sign higher 
power terms — g^x 3 — g^x A — ■ ■ ■. As explained above the cubic term corresponds to 
the gluing of triangles, the quartic term to the gluing of squares etc. The critical 
behaviour is universal as long as the coupling constants #3, #4, ... are positive and is 
uniquely determined by the quadratic nature of the maximum of V(x) for x > 0. We 



9 + Ah 2 
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(15) 
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will not consider the more general situation, where the coupling constants g$, #4, ... 
can have different signs. In that case a fine-tuning of the g^s with different signs will 
result in different critical behaviour whenever the first extremum of V(x) for x > 
will be of higher order than two. Such critical behaviour corresponds to non-unitary 
matter coupled to 2d-gravity. The non-unitarity is not unexpected in view of the 
surface representation, since it corresponds to gluing certain polygons to the surface 
with negative weight. We will here restrict ourselves to the simplest potential, 
which also has the nicest interpretation in terms of geometry: the construction of 
triangulated surfaces: 

V(x) = -x 2 - g 3 x 3 

It is convenient by a simple translation and rescaling to use it in the form 

3 

v(x) =gx~Y (20) 

Exploiting the assumed U(N) invariance of the vacuum we follow the classical treat- 
ment in [TJ| and diagonalize the hamiltonian (|1|) with respect to the eigenvalues 



Xi, i = 1, ..,N of the rescaled matrix <p = (f)/\/ N 



H fAA = -TZTZ + VfpW ( 22 ) 



N 

H FP [<t>] = N^Hfplxi] (21) 

i=l 

PPdx 2 

V fp (x) = ±(g-x 2 ) 2 + x (23) 

For the simple potential (|20|) we see that H FP [(j)} just becomes the sum of N non- 
interacting single particle hamiltonians. The particles behave as fermions since the 



expectation value of any U (N) invariant observable , calculated according to ( [3~8| ) 



is 

N 

(QM) = / Il^ A2 ({^})^o({^})Q({^}) (24) 
J i=i 

where A({xj}) comes from the integration over the angular part of the variable, 
and is given by the Vandermonde determinant: 

A({x l }) = H(x i -x j ) (25) 

i<j 



In this way, as already observed in |I8|, the function 



^({x,})^ A({xi})* (M) (26) 
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becomes the totally antisymmetric groundstate wave function of a fermionic system 
of N non-interacting particles and can therefore be written as the Slater determinant 
of the N lowest single particle eigenfunctions of Hf p [x\: 



1 



(27) 



where 



(28) 



and where the ground state energy of the total system is 




N(Ex + --- + E N ) 



(29) 



In the context of the 1/iV-expansion the total energy Eq is zero. As explained 
above the ground state energy of Hpp is trivially zero when the original problem is 
well defined and fljjD is normalizablef]. Even if the original matrix problem looks 
unbounded (it involves 3 terms, wrong sign 4 terms etc) it is actually well defined 
to any finite order in the 1/N expansion as long as the coupling constants are less 
than the critical coupling constants which determine the radius of convergence of 
the sum of diagrams to leading order in the 1/N expansion. This is reflected in the 
fact that the distribution of eigenvalues of the matrices to the leading order in N 
is restricted to a finite interval. The wave function, as given by fll4|) , is therefore 
normalizable in this leading approximation, and in fact to all finite orders in 1/N. 
Using the stochastic regularization one is able to go beyond the 1/N expansion and 
in addition one is able to extend the theory beyond the critical value of the coupling 
constant. When the coupling constant is beyond the radius of convergence the 
matrix model is not defined even in the planar approximation. But again stochastic 
regularization will provide us with a definition of the theory in this region of coupling 
constant space. 



2 This is a general property of semi-definite operators of the form (|l6|), as is well known from 
supersymmetric quantum mechanics. However, we feel that supersymmetry in this context is 
somewhat of a red heering. First, the supersymmetry alluded to in the context of matrix models 
has not yet been given a useful physical interpretation, and next the phenomenon is more general. 
As an example we can mention the case of ordinary 4d Yang-Mills theory where the operator R(x) in 
temporal gauge is iS/SA(x) — iB(x), and an exact solution to the functional Schrodinger operator 
corresponding to E = is exp(—N cs [A]), where iV cs [A] denotes the Chern-Simons number of 
the gauge configuration. However, this solution is not normalizable, as is the case for the true 
groundstate which has Eq > 0. 
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3 The large-N limit 



Since the system of TV non-interacting particles is fermionic and since we have to 
perform the summation over the first N eigenvalues, let us denote the energy of the 
iV'th level by Ep (the Fermi energy). From (^) we see that large N corresponds to 
the semiclassical limit h = 1/N 0. We expect therefore the WKB approximation 
to be good for all but the lowest eigenvalues, unless special circumstances occuiQ. 
According to the Bohr- Sommerf eld quantization rule the phase space of the classical 
theory is related to the n'th energy level of Hf p by: 

27rnh = J dpdx 6[E n - (p 2 + V fp (x))} (30) 
where h = 1/N. This gives the following relation: 

rx r (E n ) J 

/ dx JE n - V fp {x) (31) 



77, 1 fX r (E n ) 

N = n 



where x r (E n ) and xi(E n ) are the classical turning points for a particle with energy 
E n moving in the potential Vf p (x). The Fermi energy is then determined by 



1 rxr(E F ) 



1 = - dx \ Ep - VfJx) (32) 

7T Jxi(E F ) V 



while the energy density is 



dn N fME) 1 
p(E) = — = — / dx . 33 
^ ; dE 27rJ Xl (E) J E - V fp (x) 



We can now simplify the expression (p24|) if we restrict ourselves to observables 
of the form 

1 1 N 

^Tr/(0) = -E/(^)- (34) 

By expanding the Slater determinant and using the orthogonality of the single par- 
ticle wave functions \n we get 

(IlV /(«/»)) = Jdxu(x)f(x) (35) 

I N roo 

u ( x ) = j^J2x 2 i(x) ; I dxu(x) = l. (36) 



i=l 



— oo 



3 We shall discuss these special regions of coupling constant space later. 
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The density u(x) can be calculated in the WKB approximation, since the eigen- 
functions Xn{%) in this approximation in the classically allowed region where the 
energy E n > Vf p (x) are given by 

*» w *> -v,:lw« cos ( N F W* 3 ^)) (37) 

The normalization is fixed by 

POO 

dx x 2 {x) = 1 (38) 



where the integration is effectively cut off at the classical classical turning points 
xi(E n ) and x r (E n ). Due to Riemann's lemma the cos 2 is replaced by 1/2 in the large 
N limit and 

1 1 rx r (E n ) dx 71 

= 2 L l{En) l En _ Vfp (x) = N p[En) m 



We conclude that u(x) in the large A" limit is 



1 N 
JV i=i 



1 fE F 1 1 , 

— / dE . = = -JE F - VfJx) (40) 

2nJv fp(x ) J E - V fp (x) ^ fP 



It is clear from ( |33D that u(x) should be given the interpretation of density of 
eigenvalues in the original matrix model. Formulae (p5|)-(|40|) give the corresponding 
derivation in the context of stochastic regularization and we see that the finite range 
of eigenvalues in the large iV limit is determined by the classical turning points 
corresponding to the Fermi energy Ep- 

Although we know from the general theorems that u(x) has to agree with the 
eigenvalue density in the original 0-dimensional matrix model, it is interesting to 
see how this comes aboutf]. From the classical work of Brezin et al. we know 
that the semiclassical eigenvalue density of a matrix model described by a potential 
V(x) is given by 



<x) = y^Y^(^) - \(V'(x)) 2 + ft-sG*) (41) 



if V(x) is a polynomial of order k. In fl41|) the polynomial P^-^x) of order k — 3 is 
fixed by requiring the support of u(x) to be connected, In the case where k = 3 we 
get 

1 



u(x) = -yJPo-V fp (x) 

71 v 



J This argument was first presented in |io| . 
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and the constant Pq is fixed by the normalization condition for u(x). A glance on 
(|32]) allows us to identify Pq with the Fermi energy in the potential Vf p (x). 




Fig. 2 The figure to the left (fig. 2a) has g > g c . For the other figure (fig.2b) g < g c . 
The dashed lines represent the Fermi energy in the two cases. 

If we return for a moment to the original matrix model given by the potential 
V(x) = gx — x 3 /3 the sum over planar diagrams which constitutes the large N 
limit has a finite radius of convergence as a power series in the coupling constant 
(l/g) 3//4 , as explained in the introduction^. This radius is determined from the 
eigenvalue distribution u(x). The power series is convergent as long as the function 



Vfp(z) — Ep, z complex, has a single cut on the real axis |18fl . This is only possible 
if the fourth order polynomial Vf p (x) has the form shown in fig. 2a and if Ep 
coincides with the second minimum of the potential, since in this case 

V fp (z) -E F = ^{z - x + f{z - x r ){z - Xl ) (42) 

where x + denotes the position of the local minimum at the far right, while xi and 
x r are the classical turning points in Vf p for the energy Ep. It is a remarkable fact 
that Ep coincides with the local minimum at x + ||. From fig. 2a it is clear that there 
will be a critical value g c below which ( |42|) cannot be realized and where x + will in 

5 It is a power series in (l/g) 3 / 4 rather than g due to the rescaling and translation x 2 — gx 3 — > 
gx — x 3 /3. 
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fact split in two complex conjugate zeroes. For this value of g c x + and x r merge to 
a single point x c where 

V^x c ,g c ) = 0, V: x (x c ,g c )=0. (43) 

We find from (1431) 



9c =2^3' Xc = \/f- (44) 
The planar expansion will be convergent for g > g c and the eigenvalue distribution 
for a given g > g c will be located between xi and x r . This is true to any finite order 
in the 1/N expansion. Note however that while the original O-dimensional matrix 
model made no sense below g c the stochastically regularized model is perfectly well 
defined even in this range of coupling constants in agreement with general results 



4 The double scaling limit 

Let us consider the situation where g > g c (fig. 2a) and show how the double scaling 
limit as a function of only one parameter, 

n 2 = (45) 

N 2 {g-g c fl 2 (45j 
appears in a trivial way as the WKB expansion of the one-particle hamiltonian Hf p . 
If we introduce a scaled variable 



z = \_^ X _^L ( 46 ) 
V 9 ~ 9c x c 



the Fokker-Planck hamiltonian (E2l) can be written as 



Ja-a ) 3 ^ 2 

H fP = V fp (0; g-g c )+ x ^ 9 ™ h fp (z) (47) 

9c 



d 2 



h f P ( z ) = -h ~n + v fp( z 'i Vg - 9c) (4* 



dz 2 

Vf P (z; Vg - 9c) = -3z + z 3 + 



3 , Vg - gc 



--Z + -Z 

2 4 



= v ds (z) + o(^g-g c ) (49) 

The constant V) p (0; g — g c ) is a second order polynomial in g — g c - The double scaling 
limit is defined as N — > oo, g —>■ g c , h fixed, and we see that the physics in this limit 
is determined by the hamiltonian 

d 2 

h ds (z) = -h 2 —+v ds (z) (50) 
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To be precise we have dropped the trivial analytic behaviour in g — g c present in (fT7p 
in the term V/ p (0; g — g c ). Likewise one has to extract a power of g — g c in order to 
get a non-trivial result in the limit g — ► g c , as is well known from the conventional 
analysis of the O-dimensional matrix models. 

The hamiltonian (|5UD admits a WKB expansion in Ti 2 . The fact that O-dimensional 
matrix models allow for such an expansion was a non-trivial result in the original 
formulation (as well as in the continuum approach using Liouville theory), but by 
the method of stochastic regularization we get it almost for free, as was also the 
case for the eigenvalue distribution and the critical point in the large N limit. 

Since the expansion in the double scaling limit is identical to the WKB ex- 
pansion of ( |50"D it is worth while to formulate explicitly the exact WKB equation 
corresponding to flSOl). The leading order was already discussed above, but we can 
write in general: 

X (z) = A(z)e ±iS{z)/h (51) 
A(z) = A^{z)+h 2 A^{z) + --- (52) 
S(z) = S (0 \z) + h 2 S w (z) + --- (53) 

and the Schrodinger equation 

h ds (z)xe(z) = exe(z) (54) 
is equivalent to the exact WKB equations: 

S' = c/A 2 (55) 



- h 2 A"A' + c 2 = (e - v ds {z))A\ (56) 

The constant c in ( (55|) is determined by the normalization of \ e . If we introduce 
R = A 2 ([56]) can be written as 

- h 2 (R':(z)R e (z) - lK(z) 2 ) + 2c 2 = 2(e - v ds (z))R 2 e (z) (57) 
or differentiating in order to get rid of c, if wanted: 

n 2 



+ v' ds {z)R e {z) = 2(e - v ds (z))R' e (z). (5f 



This is the so-called non-perturbative equation first derived in JTJJ by means of 
the Dikii-Gelfand equation. We see that it is nothing but the WKB equation for 
A 2 {z) and in our opinion it seems to have no advantage compared to the original 
Schrodinger equation (|54|). 
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At first sight it appears as if we have gained nothing if we take the double 
scaling limit as in (|50|). The problem of dealing with the original unbounded action 
has been replaced by the problem of how to deal with the hamiltonian (0), which is 
unbounded from below. However, hamiltonians unbounded from below are perfectly 
respectable objects from a mathematical point of view, as we will review in the next 
section. They have not played a significant role in quantum mechanics, due to the 
lack of situations where they appear in a natural way. It is somewhat paradoxical 
that one has to go all the way to quantum gravity to find such a situation. 



5 Hamiltonians unbounded from below 



Consider a hamiltonian like (50), where the potential is unbounded from below. To 
be more precise we will write 

h = -£-2+<y) ( 59 ) 

where we assume that v(y) — ► — oo for y —>■ — oo so fast that 

dy 



f 

J — cx 



v(y)\ 



< oo (60) 



This condition means that a classical particle, once it is not trapped in local energy 
minima of the potential v(y), will move to y = — oo in a finite time since a classical 
particle with energy e and hamiltonian p 2 + v(y)wil\ move from y t to — oo in time 

T = f - T ^= (61) 

-v(y) 



provided e — v(y) > for all y < y t . The dynamical problem is not classical 
complete, in the sense that we need to specify some boundary conditions at y = — oo 
if we want to be able to address dynamical questions ranging over all times. From 
this point of view the mathematical aspects of the situation are not much different 
from the considerations of dynamics in a finite box. The situation transfers to the 
quantum mechanical well. Had fl59j) been defined in a finite box we would 

have been forced to impose a boundary condition for each wall in the box in order 
to get a self-adjoint operator. For each wall there would be a one-parameter family 
of self-adjoint extensions of the symmetric operator defined by (]59|) and acting on 
functions with support which does not include the coordinate of the wall itself. The 
spectral theory of such operators is well known. In the case of two walls and no 
singularities of v (y) the spectrum is purely discrete and extends to infinity for any 
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of the self-adjoint extensions. What is less known to physicists is that this result 
extends to potentials unbounded from below and satisfying (|60|). In fact the usual 
situation of walls with imposed boundary conditions can be treated as a special case 
of the potentials unbounded from below, and although this might seem somewhat 
perverse from the point of view of physics, it is not unnatural from a mathematical 
point of view, where both situations can be classified as Weil's circle-limit case of the 
Sturm-Liouville theory of second order differential operators (see for instance [l9j for 



a recent discussion). The important conclusion is that the spectra of the self-adjoint 
extensions of ([59D are purely discrete and extend from — oo to +00, provided is 
satisfied and the potential behaves in a similar way for y — > +00 or v (y) — > 00 for 
y -> 00. 

It is not the purpose here to describe in any detail the Sturm-Liouville theory 
for potentials unbounded from below, but the important points are easily explained. 
Suppose we want to solve the differential equation 

hip = eip (62) 



where h is given by (|59f) and ( p0|) . Far to the left the potential goes to —00 and the 
WKB-approximation becomes an excellent approximation. In this region we can 
write: 

1 [ ( rvt(e) 



ip(y; e, a) ~ ^ _ ^^1/4 cos j yj V e ~ v ( x ) dx + a j j ( 63 ) 

where a is an arbitrary angle and yt{e) denotes the classical turning point for a 
particle coming from the left (we assume that e is so small that this point exists). 
We see that any solution to ( |62"D is square- integrable at —00 thanks to the condition 
(|6C|). However, this class of functions is too large to constitute the domain of a 
self-adjoint version of h since a partial integration leads to a term: 

# 2 (z/;e 2 ,a 2 ) #*(|/;ei,Q!i) 

V>iO/;ei,ai) ip 2 {y; e 2 , a 2 ) (64) 

dy dy 

which does not go to zero for y — > —00. Special choices of a ensure the convergence, 
namely 



rjft(e) 

a(e) — — I dy 



—00 



Je-v(y) - Je -v{y) +/ Je -v(y) (65) 

J Jyt[e) 



where eo > e (and where we again assume that eo is chosen such that yt(eo) exists). 
We see that the integrals in (|65| ) give a well defined way to write the formal difference 
in WKB phases corresponding to energies e and eo: 



me 



yt(e ) 1 ryt(e) 



dy\Je -v(y)- dy^Je-v(y) (66) 
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With this choice all functions in (p3f) have the asymptotic behaviour 



and for fixed eo ( p4|) goes to zero for y — > — oo. The different self-adjoint extensions 
will be characterized by different choices of the parameter eo- Due to the asymptotic 
behaviour ( |67D of the wave functions it is possible to characterize the unbounded 
hamiltonian as a sequence of ordinary bounded hamiltonians [^(J . For — y sufficiently 



large all wave functions VKz/i e > e o) will vanish at points y n given by: 



y/-v{y r . 



Consequently one will get the same result if one replaces the unbounded potential by 
a potential v n (y) cut off by an infinitely high wall at y n , provided |e| is not too large. 
In this way the self-adjoint hamiltonian, unbounded from below and characterized 
by the parameter eo, is reached by a sequence of ordinary hamiltonians, bounded 
from below, and corresponding to the potentials v n (y). 

This discussion, based on the WKB approximation, can be made mathematical 
rigorous, but it is worth emphasizing that it can be formulated independently of the 
WKB expansion. Let e be given and let i[>\(y;e) and ^2(2/; e) be two independent 
solutions to (|62D with a wronskian (which of course is independent of y) 

w = ip' 2 (y)ipi(y) - ip 2 (y)ip[(y)- (69) 
If we define a class of functions by 

f(y) = if>i(v',e) (ci + — [ dx tfj 2 {x;e)g(x)\ + 

I W J -00 J 

^2(2/; e) (c 2 - — / dxif> 1 (x',e)g(x)\ (70) 

where g is any square-integrable function, it is readily seen that 

hf = ef + g and f(y) -> Ctfj)\(y\ e) + 02^2(2/; e) for y -»• -00. (71) 

Since / is square-integrable, h maps / into a square-integrable function, and using 
the asymptotic behaviour ([T!]) in ( p^ ) it is seen that the class of functions defined 
by ([7t]) constitutes the domain of a self-adjoint version of h provided the socalled 
(complex) limit numbers c\ and c 2 satisfy 

ci cos «o + C2 sin a = (72) 
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for some real olq. In this way we recover the WKB results above. 

If we assume that the potential also satisfies a condition like (j60|) for y — > oo we 



can repeat the discussion above and get a relation like (72), just with other constants 
c[, d 2 and a' Q . In order to solve the eigenvalue equation we have to choose g equal to 
zero in (|70|) and we get a matching condition at (say) y = for / and /'. For fixed 
«o and a' the matching can only be satisfied for certain discrete values of e, leading 
to a quantization of the energy eigenvalues in much the same way as for a particle 
in a box of finite width. If the potential instead goes to +00 for y — ► 00 we also 
get a purely discrete eigenvalue spectrum. The solutions to (|62"D will in this case 
consist of exponentially growing and exponentially decaying parts for y —>■ 00. The 
linear combination of ipi(y;e) and ^2(2/; e) dictated by ( [72] ) will in general contain 
an exponentially growing part. Only for special discrete values of e will the wave 
function be exponentially decaying and therefore square- integrable for y — > 00. This 
last situation will be the one of interest to us. 

6 Application to Vd s (z) = —3z + z 3 

Let us apply the formalism of unbounded hamiltonians to the hamiltonian ( pTp and 
compare the results with similar results obtained by using the full Fokker-Planck 
potential. We first discuss the quantization of energies. The potential has a local 
maximum at z = — 1 and a local minimum at z = +1 with a value Va s (z = +1) = —2. 
The same value of the potential is obtained if z = —2, which therefore is the classical 
turning point when the energy coincides with the local minimum, as is the case for 
the Fermi energy in the WKB limit for the full Fokker-Planck potential Vf p (z) given 



in (fEJ). Let us denote the Fermi energy of the full Fokker-Planck hamiltonian (]48 
by ej. It is trivially related to the E F defined previously by 

(a- a W 2 

E F = Vfr(0; g - g e ) + x\ e f . (73) 

As mentioned above the spectrum for hd s {z) is purely discrete. If h is small and 
e < ej we know that the solutions to the Schrodinger equation are well approxi- 
mated by the WKB solutions ( |63| ) to the left of the classical turning point z t (e). To 
the right of the turning point it is given by similar expressions which fall off or in- 
crease exponentially. The WKB matching condition to the exponentially decreasing 
solution is 

^; e) ~ {e _ v 1 ds{z))1/i ™{lr (e) ^-v ds (z') dz> - tt/4} (74) 
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and the quantization of energies in the WKB expansion comes about by comparing 
with the phase requirement fl6"5|) which was needed in order that we had a self-adjoint 
hamiltonian. 

If we let e' , the parameter characterizing the self-adjoint extension in the WKB 
language (not an eigenvalue), be of order one, we get for the n'th eigenvalue below 

rz t {e n ) / J , \ rzt(e' ) , 

-vr^(n + 3/4) = / dz [Je n - v ds (z) - Je' - v ds (z) - / dzJe' -v ds (z). 

J -co \ J J z t (e n ) 

(75) 

We can solve this equation in two regions (assuming as usual that h is small). 
Assume first that e' — e n is of the order of ft, and denote by e the first eigenvalue 
below e' . We get: 

e n = e - 2hun + o ((ftn) 3/2 ) (76) 

where 

rz t (e ) | 
U = 7l/ dz = (77) 

J -°° \Je -v ds (z) 
has the interpretation as the cyclic frequency in the classical motion out to infinity 
in a finite time (see (|6~ID). The only arbitrariness present at the energy levels close 
to the Fermi energy is therefore an all-over displacement of the energy levels. If we 
choose e' such that eo = e/, the Fermi energy, we get 

uj = u f = V3, (78) 

which is in agreement with the similar calculation done in the full Fokker-Planck 
potential: 

e n = e f - 2hu f n + o {{hnf' 2 ) +o((g- g c ) 1/4 ) . (79) 

The main source to the difference between flTBD and (|79|) is the difference in escape 
time to infinity in the potential v ds (z) compared to the time it takes to reach the 
left turning point in the Fokker-Planck potential Vf p (z). This difference is of the 
order (g — g c Y^ and vanishes in the double scaling limit. 

The above results are even more manifest if we calculate the energy density for 
small h where we expect (|73p to be a good approximation and where the eigenvalues 
are dense 

dn ds 1 r Zt( - e > dz , , 

1 (80) 



d\e\ lixriJ-co yje-v^z) 

We see that all reference to e' Q drops out. Let us compare this result with the corre- 
sponding result for the full potential Vf p (y). From ( |33D we get by simple rescaling: 



drif p 1 r 2 '-^) dz 



d\e\ 2nh J tl (e) J e - Vfp (z) 
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where z r (e) and zi(e) denote the right and left turning points. We have z r (e) = 
Zt(e) + o(yfg — g c ) and one can check that 

= + o((g - g e )W) (82) 

provided e is not too close to the bottom oivf p (z), i.e. provided e >> —(g — g c )~ 3 / 2 . 
Again we see that the difference between the eigenvalue densities of ha s and hf p 
vanishes in the double scaling limit as (g — gc) 1 ^, and that this is valid for a much 
larger energy range than indicated by ( |75| ) and flT5|). 

The other region where we can solve (|75|) is for — e n » — e/ ~ 2. In this region 
we get 



where the constant c is 



We can write (|83|) as 



c I 



(l+o((nh)-^)) (83) 



-y 3 -V-i-y 3 ) + y^-?/ 3 ( 84 ) 



,6/5 (^-^ c ) 3 / 2 1 J 



and even in the case of the the Fokker-Planck potential will this expression for the 
energy be valid all the way down to the iV'th level below the Fermi energy. However, 
the sum of these large negative eigenvalues does not contribute to the non-trivial 
scaling since we have, using (|85|) 



N 



AE~N(g-g c )V 2 Y,en~N 2 . (86) 



Although the energy eigenvalues of the self-adjoint extensions agree with the eigen- 
values of the full Fokker Planck potential in the double scaling limit up to an all-over 
displacement we need reference to the full Fokker Planck potential in two ways if 
we want to calculate physical observables like (Tr <p) or more generally (p4|) . We 
need it to define the top level (i.e. the Fermi energy) and the bottom level in the 
summation (pop, but we also need it to define a cut-off for large negative z when we 
calculate expectation values of operators like (Tr <fi). The reason is that although 
the density u(z) defined from ( |3"6"D is integrable out to z = — oo Uf p (z) and Ud s (z) 
will differ when we go beyond the left turning point of for the full Fokker Planck 
potential. Beyond the left turning point, 

*< e > ^ -70S (87) 
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Uf p (z) falls of exponentially. It even vanish beyond z\{e) in the large N limit. This 
is not the case for Ud s (z) and integrals like f A z n Ud s (z) diverge for A — > — oo. If we 
introduce a cut-off 

A — (88) 

\j9~9c 

we get however the correct scaling behaviour of our observables. The only features 
we need to borrow from the full Fokker Planck potential in order to get the correct 
scaling behaviour of observables in the double scaling limit is thus a Fermi energy 
of order one below the local maximum of Vd s (z) and a cut-off A in the negative 
z as given by fl88|). Note that such a cut-off by an infinite wall is by no means 
unnatural for the self-adjoint extensions since their wave functions converge to the 
same oscillating function (|6?|), independent of the eigenvalues e„, for z — > — oo. We 



can now put the wall at one of the zeros of the oscillating function ( 67[) and this 
will not influence the spectrum or eigenf unctions. If we are only interested in the 
non-trivial scaling behaviour of the physical observables which is associated with the 
asymptotic WKB expansion in h 2 all of the self-adoint extensions can be used. If we 
however are interested in the full non-perturbative contributions which survive in 
the double scaling limit, but which cannot be expressed as powers of h 2 , the different 
self-adjoint extensions will differ, and we have to use either the full Fokker Planck 
potential or the unique self-adjoint extension which has an energy eigenvalue exactly 
equal to the Fermi energy e/ of the full Fokker Planck potential. This will be clear 
in a moment when we calculate the lowest non-perturbative correction in the limit 
of small h. 

Let us finally address the non-perturbative corrections which are not to be be found 
as powers of h within a systematic WKB expansion. When g > g c we have the 
situation in fig. 2a. We have two wells. A large one (unbounded from below in the 
double scaling limit) and a smaller one above the Fermi energy to the right. In the 
limit of small h the tunneling through the barrier will be exponentially suppressed 
as e~ r ^ h , where T will be calculated below, and to a first approximation we have 
independent states in the left well and the right well. However, the situation is 
remarkably similar to that of a symmetric well, even if there at first sight is nothing 
symmetric about the situation. First we note that the level spacing is the same on 
both sides of the barrier. From ( |T?D and (|78| ) we get that the level spacing in the 
left well is given by 2hujf where Uf — \^3. The level spacing in the right well is 
determined by expanding hd s (z) around the local minimum at z — 1 and we get 
with z — z — 1: 

h ds (z) = -h 2 -^ + 3z 2 + z 3 (89) 
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and the harmonic frequency which determines the lowest lying levels in the semi- 
classical approximation is just uj = y/3. Not only is the level spacing the same on 
both sides of the barrier, but in addition the Fermi energy in the well to the left 
coincides precisely with the lowest energy in the well to the right. 



Fig. 3 The energy levels in the two wells. The lowest energy in well to the right 
coincides with the Fermi energy in the well to the left if this energy is calculated 
according to the WKB prescription. 

Until now we have determined the Fermi energy by ([^), which in scaled variables 
reads 

UN = — dz Je f - Vf p (z). (90) 

IT Jzi(e f ) 



As remarked above (see (|42|)) this implies that e/ = —2, i.e. e/ coincides with the 
local minimum of the right well. However, to be precise one should in the context 
of WKB expansion replace N by N + | and this leads to a shift of the Fermi energy 
such that it agrees with the lowest energy in the right well. We have now a situation 
identical to the standard double well. By tunneling the degenerate energy levels 
split in two and the lowest one will now have to be identified with the Fermi energy. 
The shift is non-perturbative in h and is given by the standard WKB formula: 

C/ - C/ - ^ . e-W (91) 

7T 

where the tunneling amplitude V is given by 

dz \jv fp (z) - e~ f « f dz \Jv ds (z) + 2 = ^r^. (92) 
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In zi is the position of the left turning point for a classical particle with energy 
e/ in the right well, while z r as usual denotes the right turning point for a classical 
motion with energy Cf in the left well. Again we see that it is possible in this 
semiclassical calculation to replace Vf p by v^s, in which case the integral can be 
performed since Vd s (z) + 2 = (1 — z) 2 {2 + z) 

The lower levels e n , n > will receive contributions which are exponentially 
small in 1/h compared to the one received by ef since they are located below the local 
minimum of the well to the right. We note that the contribution (j9ip is identical 
to the "non-perturbative" ambiguity which comes from the Painleve equation of 
ordinary gravity (but in the stochastic approach it is of course not an ambiguity). 
This observation was first made in refs. |^T], |22| and is in agreement with the fact 
that the WKB expansion of the stochastically regularized model agrees with the 
genus expansion of the original matrix model of pure 2d gravity. 



7 The strong coupling expansion 

The theory outlined above offers the possibility of a convergent strong coupling 
expansion. It might be useful to recall a similar situation in the case of the ordinary 
anharmonic oscillator in quantum mechanics. Let the hamiltonian be given by 

d 2 

H = -— + x 2 + gx\ g>0. (93) 
dx A 

The ordinary perturbation expansion around the harmonic oscillator is only an 
asymptotic expansion. The reason is that g x 4 is not a small perturbation for any 
value of g > and clearly g — > —g changes drastically the nature of H. For this 
reason the expansion in powers of g is not a convergent power expansion, but only 
an asymptotic expansion. It is, however, possible to analyze the strong coupling 
region g — > oo by a simple scaling argument: 

H(x) = \- 1 / 2 H(y), y = g 1/6 x, X = 4n ( 94 ) 



g2/3 



H(y) = -f^ + y 4 + \y 2 (95) 

From this we deduce that H(y) has a strong coupling expansion, which is analytic 
in A = l/g 2 / 3 . In fact the potential V(y) = y 2 satisfies the standard requirement for 
being an analytic perturbation of 
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namely that V is i^-bounded, i.e. (1) the domain D(V) D D(H ) and (2) ||V^|| 2 < 
a\ \H ip\ | 2 + b\ \1pW2 for some a and b and all ip G D(H ). For an eigenvalue E n of the 
original hamiltonian H we can therefore write: 



A- 1/2 £ c nk \ k 



(96) 



k=0 



where the power series has a finite radius of convergence. 



The situation is very similar in the case of the hamiltonian hf p (z). In sect. 4 we 
introduced a scaling which was designed to make contact with the asymptotic genus 
expansion of ordinary 2d quantum gravity, i.e. the limit H — > 0. But we can perform 
the same scaling argument as just outlined for the anharmonic oscillator and derive 
an expansion for large h 2 , i.e. in the strong coupling regime. If we introduce scaled 
variables 

1 9 ~9c 



5 



V: 



we can write 



A 9c 

h fp (z) = \- 3/ % P (y) 



N\ 4/5 
2 ) 



(97) 
(98) 



hf P (y) = 



d 2 

■— + v fp (y,X,N) 



(99) 



/ 2 \ 2/5 



= v ds (y, A) + ( — 



2/5 



3A 2 1 4 

- Y y + - A y 

3A 2 1 4 

Y y +- A y 



(100) 



and the double scaling limit is obtained as before for N — > 00, g — > p c , but A fixed. 

The non-trivial information is contained in the Hamilton function hf p (y) which 
in the double scaling limit goes to 



hds{y) = 



d 2 



+ Vds 



d 2 



+ y 3 - 3Ay 



(101) 



dy 2 dy 2 

and we have formally the same situation as for the anharmonic oscillator: — 3Xy 
looks like a small perturbation with respect to y 3 . It is thus natural to expect a 
strong coupling expansion of the energy eigenvalues of hd s (z) of the form: 



e n (h) 



A" 3/2 £ c nk \ k 



k=0 



k=0 



(102) 



where the series is convergent. The domain of the perturbation v(y) = —3\y does 
not include the domain of ho(y) = —d 2 /dy 2 + y 3 so we have no rigorous proof of this 
conjecture, but one should keep in mind that the requirements (1) and (2) mentioned 
above are only sufficient conditions, not necessary conditions, for analyticity. 
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8 Discussion and conclusion 



We have shown that many of the results of the matrix models are easily and trans- 
parently derived by means of stochastic quantization. The asymptotic expansion, 
called the double scaling limit, is nothing but the WKB expansion of the Fokker- 
Planck hamiltonian. In addition stochastic quantization provides a non-perturbative 
definition of 2d quantum gravity. To be entirely correct one should first calculate 
expectation values of observables using the full Fokker-Planck hamiltonian and af- 
terwards take the double scaling limit N — > oo, g — ► g c , with ft ~ N~ l (g — g c )~ 5 ^ 
fixed. If we reverse the procedure the double scaling limit of the Fokker-Planck 
hamiltonian results in a hamiltonian hd s which is unbounded from below and which 
has a one-parameter family of self-adjoint extensions. One member is picked out 
by the requirement that the Fermi energy should coincide with the corresponding 
energy of the full Fokker-Planck potential. We mentioned the possibility of a strong 
coupling expansion, similar to the one of the anharmonic oscillator. 

It is interesting to compare the situation described above with the similar situ- 
ation present for 2d gravity coupled to a scalar field (i.e. c = 1). In this case we 
again have a matrix model description which, in the double scaling limit, results in 
a Schrodinger eigenvalue equation: 



d 2 

j-2 u 2 



ip n = e n ip n . (103) 



where h' 1 = N(g c - g) is fixed ||, H, [2S|. A gain the ground state is fermionic 
and given by the Slater determinant of the N eigenfunctions counted from the top 
of the upside-down quadratic potential in ( |103j ). One could be tempted to apply 
the methods outlined above to this hamiltonian, which is unbounded from below, 
but condition ( j60| ) is not satisfied for this potential, and the problem is classical 
complete: It takes an infinite time for a classical particle to escape to infinity. In 
the terminology of Weil it is a point-limit case of the Sturm-Liouville theory of 
second order differential equations. The point spectrum is empty and there is a 
unique self-adjoint extension. If we cut off the potential at ±L we get a sequence 
of Hamiltonians which are bounded from below and which convergences to the self- 
adjoint extension of the lhs of ( 103j) for L — > oo. The (generalized) eigenfunctions 



are parabolic cylinder functions and they allow for the calculation of both local |26 



and non-local [^] observables within the WKB approximation, but although the 
self-adjoint extension of hamiltonian given by ( |103| ) is unique and defined beyond 
the WKB expansion, it is still not clear whether it allows a strong coupling expansion 



for small h as conjectured in |24 
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